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We calculate the low energy elementary excitations of a Bose-Einstein Condensate in an effective 
magnetic field. The field is created by the interplay between light beams carrying orbital angular 
momentum and the trapped atoms We examine the role of the homogeneous magnetic field, 
familiar from studies of rotating condensates, and also investigate spectra for vector potentials with 
a more general radial dependence. We discuss the instabilities which arise and how these may be 
manifested. 
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I. INTRODUCTION 

Quantum degenerate gases are in many ways the ideal 
quantum many-body system. In an experimental situa- 
tion they afford an unprecedented level of control over 
the system parameters, such as the strength (and even 
sign) of the interaction between the atoms, the geometry 
of the external trap and the properties of the lattice in 
which the atoms are loaded. It is no surprise therefore 
that Bose-Einstein Condensates (BECs) and degenerate 
Fermi gases are often used as a laboratory to study a host 
of phenomena from many different areas of physics. This 
is especially true in condensed matter physics; for exam- 
ple, ultracold atoms in an optical lattice can be stud- 
ied using the Hubbard model Similarly, a system of 
trapped fermions tightly confined in one direction invites 
obvious comparisons with the 2D electron gas [3|. 

Without doubt, some of the most striking effects in 
solid state physics are observed when an external mag- 
netic field is applied to a collection of charged particles. 
Well known examples include the quantum Hall effects 
in 2D electron gases and the Meissner effect in Type II 
superconductors. As the atoms forming quantum gases 
are electrically neutral, it is not obvious at a first glance 
how they might be used to study such effects. 

The solution lies in the ability to create artificial 
magnetic fields. For example, rotating the system and 
studying it in the rotating frame is analogous to study- 
ing charged particles in a homogeneous magnetic field 
0) IE Hi- Alternatively, lasers can be used to alter the 
state-dependent tunneling amplitudes of atoms in an op- 
tical lattice to simulate an effective magnetic flux @, H|- 

A recent proposal involves the adiabatic motion of 
lambda-type three level atoms interacting with laser- 
fields which create a non-degenerate dark state, that is 
an eigenstate of the atom-laser interaction. It has been 
shown that if the atoms interact with a pair of laser 



beams possessing a relative orbital angular momentum 
[l(| , then an effective vector potential appears in the 
effective equation for the atomic wavefunction P, [T]| . 
The corresponding effective magnetic field created is en- 
tirely dependent on the form of the incident light, so that 
by appropriately choosing the light's phase and intensity 
we can control both the strength and shape (homoge- 
neous or inhomogeneous) of the effective magnetic field. 
The inherent flexibility of the system allows for wide- 
ranging studies into the magnetic properties of both de- 
generate Bose and Fermi gases, and could provide insight 
into gauge theories in general. 

It is therefore pertinent to gain an understanding of 
how the fundamental properties of the gas may be mod- 
ified in the presence of artificial magnetic fields. A com- 
plete analysis must include the excitations, which deter- 
mine the dynamical behaviour of the system under weak 
perturbation and are crucial in determining its superfluid 
properties. Of particular interest are the lowest energy 
(or elementary) excitations, which are collective in na- 
ture. In this paper we calculate the spectra for a trapped 
2D BEC in both homogenous and non-homogeneous 
magnetic fields, which are created as described in 
Two-dimensional quantum gases have recently attracted 
a considerable interest in connection with the Kosterlitz- 
Thouless transition [12j and the quantum Hall effect in 
clouds of ultracold atoms (see [3[ and references therein). 

The paper is organised as follows: in section II a brief 
description of the model is given and in section III we 
outline how the excitations are calculated. As the inter- 
action between the light and atoms introduces two ef- 
fective potentials - a vector potential and also a scalar 
potential - and both have a significant role to play in 
the dynamics and excitations, we present our results in 
two parts. In section IV we assume the external trap 
has been chosen to counteract the effect of the additional 
effective trap so that the potential felt by the atoms is 
completely harmonic. This allows us to isolate the role 
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are, respectively, the effective vector and trapping poten- 
tials. The external trapping potential for the dark state 
atoms is 



FIG. 1: Method for creating effective magnetic fields in de- 
generate atomic gases. On the left is the level scheme for 
the A-type atoms interacting with the resonant probe beam 
Q p and control beam Q c . On the right is a schematic repre- 
sentation of the experimental setup with the two light beams 
incident on the cloud of atoms. The effective magnetic field 
is generated if there is relative angular momentum between 
the beams. This will occur, for example, if the probe field is 
of the form Q, p ~ e 1 ^, where each probe photon carries an 
orbital angular momentum hi along the propagation axis z, 
and fi c is independent of the azimuthal angle. 
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where Vj is the trapping potential for the atoms in hy- 
perfine state j (j = 1, 2). The atomic interactions, which 
involve collisions between atoms both of the same and 
of different electronic state, can be described by a single 
parameter [1| 
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of the magnetic field alone on the excitations. Then in 
section V we include the full effective trapping potential 
terms and study the excitations numerically. Finally in 
section VI we discuss and summarise the main results. 



II. THE MODEL 

We consider a system of three-level atoms charecter- 
ized by two hyperfine ground levels |1) and |2) and an 
electronically excited level |3) interacting with two co- 
propagating resonant laser beams in an EIT configura- 
tion (figure 1). The probe beam, which has coupling 
strength f2 p and is allowed to have angular momentum 
£H per photon along the z-axis, drives the transition 
|1) — > 1 2), whilst the control beam has coupling strength 
f2 c and is concerned with the transition |2) — > |3). These 
absorption paths destructively interfere to suppress tran- 
sitions to level |3), driving the atoms to the dark state 

— — where 



superposition of levels |1) and |2): \D) 



u. 



3 lS and S is the relative phase between 

the probe and control beam. If the atoms in the dark 
state form a BEC, then the coupling between the light 
and the atoms introduces an effective vector potential 
into the mean-field equation for the atomic wavefunction 



ih^ = ^(zhV + A)H + V(r)y + g\y\ 2 y, (1) 



where 
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Here, gij — 47r ^ a 'j , where ay is the s-wave scattering 
length between atoms in the levels i and j = 1,2): 
a,jj is the scattering length of atoms in the same electronic 
state and a\2 = a^i corresponds to collisions between 
atoms in different electronic states. Note that in general 
the infractions can depend on position since C is position 
dependent. However, if £ is small or alternatively if the 
inter and intra-spieces scattering lengths are equal then 
g is approximately constant throughout the condensate. 
We shall assume this to be the case for the remainder of 
this paper. 

For the purposes of our analysis it is convenient rewrite 
Eq. ([!]) in the form 
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where V(r) = V + Wf • F° r the examples considered 
V • A = so this equation is equivalent to 

As can clearly be seen from Eq. ([2]), a non- vanishing 
effective magnetic field is created if there is relative or- 
bital angular momentum between the two light beams, 
such that the phase of the dimensionless ratio ( is given 
by S = £(/>, where <f> is the azimuthal angle. The shape 
of the effective vector potential is controled by the in- 
tensity ratio of the probe and control beams. We choose 
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= a ° ' w h ere the dimensionless pa- 

rameter ctQ is the ratio of probe to control beam at a 
characteristic radius r = R which is chosen to be larger 
than the radius of the trapped cloud. The exact forms of 
the resultant effective vector and trapping potential are 
then 
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In order for the adiabatic dynamics to hold, so that 
the atoms remain in the dark state for a typical BEC 
lifetime, requires, typically, the ratio |£| 2 < 1 and hence 
ao < 1 |]. If ao << 1 we can safely replace Eqs and 
flU) by the approximate potentials 



V = V P . 
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and the effective interaction strength (eq. O is constant 
throughout the cloud: g w gn. We obtain Eqs. (jSJ) and 
([ID)) simply by expanding eqs (J7J) and (jHJ) and neglecting 
all higher order terms in ao, except of course the 2nd term 
on the right hand side of (fTUj) . because this is multiplied 
by an i which can in principle be arbitrarily large. 

The freedom to choose the form of the spatially- varying 
effective magnetic field relies on the ability to shape the 
intensities and phases of the incident laser beams. Recent 
advances in light beam shaping, using for instance spa- 
tial light modulators, makes it possible to consider truly 
exotic light beams The tightly confined or two di- 
mensional gas offers in this respect a clear simplification: 
shaping a light beam in a plane is much less restrictive 
than in three dimensions, although 3D light shaping is 
certainly possible. 



III. CALCULATING THE EXCITATIONS 

We calculate the excitations of this system by consid- 
ering small time-dependent variations of the condensate 
wavefunction around the ground state \&(r), writing the 
wavefunction as 

#(r, t) = [*(r) + u(r)e- lut - v* {r)e iuJt }e~^ t/h . (11) 

Substituting this into eq. © and keeping only linear 
terms in u and v we obtain two coupled equations anal- 
ogous to the Bogoliubov - de Gennes equations [l4[ 

{£ -^i)u + g\y\ 2 (2u-v) = (n + E)u (12) 
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\)v + g\*\ 3 (2v-u) = Qi-E)v, (13) 



where Cq 
satisfies 



-^rrV 2 + Va,nd the (non- vortex) ground state 
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£ + fl|*r)* = A** 



These are the equations we solve to calculate the eigen- 
frequencies to and eigenenergies E — hu> of the excitation 
modes u and v, which satisfy the normalization condition 



dv(u a u* b - v a v* b ) = 5 ab . 



(15) 



In deriving Eqs. (|E^ - (|T4|) we have allowed the ex- 
citations to have angular momentum mh (where m = 
0, ±1,±2...) by transforming u and v such that u — > 
we"™* and v — ► ve lm ^ '. We assume that the ground state 
ip = ip(r) only; this is equivalent to saying that ip does not 
correspond to a vortical state. Note that a more riogor- 
ous approach for deriving the excitation equations would 
be to diagonalize the many body Hamiltonian in the Bo- 
goliubov approximation, expressing the fluctuation op- 
erator in terms of quasiparticle operators [HI]. The re- 
sultant equations which must be solved are identical to 

HEME). 

The excitation spectrum of eqs (fT2]) - (|T3"| can be cal- 
culated in the spirit of the Thomas-Fermi approxima- 
tion where the effective trap potential (which we assume 



to have a harmonic ('■ 



component) and the repul- 



sive mean-field interactions provide the dominant en- 
ergy scales. In this regime /i « n max g >> hu>t, where 
n-max is the maximum condensate density. We rescale 
the radial coordinates as p = £ and R = where 



L = ( — ^ is the characteristic length-scale of the 

harmonically trapped condensate with trap frequency u>t- 
Introducing the dimensionless eigenenergies e v = and 

I*! 2 

the dimensionless density n = -J — ■ — , we can re- write our 
equations as 



C Qd u + (2u - v)n = l + 2 7 (e + 7^) u (16) 

£ 0d v + (2v - u)n = [l-2 7 (e + 7 ^)] v (17) 
d^ + ml> = il>, (18) 

with £ 0d = — 7 2 V 2 + Vd where, in the Thomas-Fermi 
regime, the parameter 7 = 4^- << 1 will allow for fur- 
ther simplification of Eqs. (|16 p - (fT5|) . The subscript d 
in Vd and Ad denotes that dimensionless units are being 
used. 

Equations (fTB|) - (fT5|) are reduced to two fourth or- 
der differential equations after introducing the functions 
f± = u ± v and substituting n from Eq. (fT5|) : 
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2(e+^)/ + , (19) 



, 1-Vd . V 2 ^ 




= 2(e + ^)/_, (20) 
P 

In the Thomas-Fermi approximation, we neglect the 
kinetic energy term 7 2 V 2 ?/> in Eq. (JT5J) to obtain the 
wavefunction 



(1 - V d ), *tf > 0, 



(21) 



*f?TF 



which can readily be substituted into Eqs. (|19|) and ([20)1 . 
When considering the low energy excitations, with wave- 
functions that vary over a scale comparable with the size 
of the condensate, we must for consistency also neglect 
terms proportional to 7 2 in (fl9|) and ([20| . Then, applying 

the ansatz /± = C± (e + ^) 3 f^^j) " Q0>, 0), 

we derive in the Thomas-Fermi regime the equation for 
the excitations of a 2D condensate with effective mag- 
netic field: 



(1 - V d )V 2 Q - (1 - V d ) 
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The relation between the normalization constants C+ 
and C- can be obtained from Eqs. (fTrJ)) , (JTTJ) and 



C_ = yC. 
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(23) 
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FIG. 2: Numerical results for surface mode (n=0) excitation 
frequency as a function of the effective magnetic field strength 
aol for different vector potentials. The top left panel corre- 
sponds to the homogeneous magnetic field. The lines cor- 
respond to the five lowest energy surface excitations, where 
the m=0 case (Goldstone mode) is represented by the low- 
est frequency mode at aol = 0, and the m = 4 the highest. 
Each mode exhibits an energetic instability beyond a critical 
olqI. The mode frequencies are purely real, so the system is 
dynamically stable when the vector and trapping potentials 
are radially symmetric. A full description of the numerical 
method is given in section V. 



that V = ^Muj t r 2 1]. If the ratio of control to probe 
intensity in the transversal plane is of the form |£ 2 | ~ r 2 , 



then the exponent v in eq. @ is 1. 
potential is then 



Mao 
"R 2 



The effective vector 



(24) 



corresponding to a uniform magnetic field in the z- 
direction i.e. B = V x A = - 2 ^"° e z . In this case, the 
solution of (|22|) is of the form W = x^ P(x)e lm ^ , where 
x = p 2 and the radial function P(x) is the solution of the 
hypergeometric equation 



x(l - x)-— + ((m + 1) - (™ + 2)x) 



<PP 
dx 2 

~e 2 - jaotme - £ ) P(x) = 



dP 

dx 



(25) 



For a physically well-behaved solution, we require P{x) 
to be convergent at x — and converge as x — ► 1. This 
yields the spectrum 



IV. HARMONIC EFFECTIVE TRAPPING 
POTENTIAL 

A. Homogeneous Magnetic Field 

By a judicious choice of external trapping potential, 
a purely harmonic effective trap can be obtained, such 



1 

E n . m = hw t (2n 2 + 2n\m\ + 2n + \m\) 2 + j^a £m 

(26) 

where n is the number of radial nodes and to the angular 
momentum quantum number. This spectrum is of the 
same form as that found by Ho and Ma [l6| for the two- 
dimensional cloud, except for the shift term proportional 
to ao£m which arises due to the effective magnetic field. 
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The uniform effective magnetic field is thus shown to in- 
duce a Zeeman-like shift on the energy levels, decreasing 
the excitation energy when I and m are opposite in sign 
and increasing it when they are the same. The actual 
mode observed should be interpreted as a superposition 
of the + and - m modes, as the effective magnetic field 
induces a rotation in modes with m/0. For example, if 
we excite the m = ±1 dipole mode (which corresponds 
to a sloshing motion along one axis) and then switch on 
the effective magnetic field the mode will start to precess 
at the effective cyclotron frequency due to the additional 
force propotional to v x B, where v is the velocity. 

Let us also highlight that as a result of our choice of 
external trap, the spectrum of Eq. (|26|) matches that of a 
rotating condensate when studied in the rotating frame, 
with rotation frequency given by fl — jjfp . The stability 
properties of this system have been studied extensively 
(see e.g. [HElIl)- 

The solutions of Eq. (|25|) are the Jacobi polynomials 
p(\™\fl)(i _ 2x ) an d from Eqs. HU) and tfH]) we obtain 



H 



An + 2\m\ 



L 2 



(1-P 2 ) 



xp ™p(\™\fi)(i_2p 2 ). (27) 



Vector Potential for which the Magnetic Field 
is Zero 



If instead the ratio of the control and probe beam in- 
tensities were constant (y = 1) we would obtain the vec- 
tor potential 



-ha -e^ 
r 



(28) 



in which case B = VxA = 0so that the effective mag- 
netic field is zero throughout the cloud (in the same way 
that a velocity field proportional to 1/r around a vor- 
tex still satisfies the irrotationality criterion for Bose- 
Einstein Condensation). An approximate energy spec- 
trum can be derived by treating the effective poten- 

noting that 
The effecive 



0. 



tial as a small perturbation in eq. 
this treatment breaks down as r - 
vector potential plays the role of an additional cen- 
trifugal potential. The solution of (|22"|) is of the form 

w = x ^m 2 +i lao lrnp^ e im4>^ where x = p 2 and p^ 

is governed by a hypergeometric equation which admits 
a physical solution convergent at the origin and as x — > 1 
only if 



-e 2 — "fao£me — (n+ -)V m 2 + 4jao£me = n 2 +n, (29) 



from which we can obtain the excitation spectrum by 
solving for e. In the perturbative regime where ejaoi << 



we find 



hujt 



hut 
~2fj, 



a$£ (m + (2n + l)sgn(m)) 



>\y- 



(30) 

gives the spec- 



where e (0) = (2n 2 + 2n\m\ + 2n 
trum when there is no vector potential and sgn(m) = 
+1, —1,0 for m positive, negative or respectively. This 
is a somewhat crude approximation but, nevertheless, it 
yields an insightful result: the effective vector potential 
has significance even if the corresponding effective mag- 
netic field seen by the atoms is zero, reminiscent of the 
Aharanov-Bohm effect. This significance is manifested 
in a shift in excitation energy levels for all modes with 
angular momentum (m 7^ 0), and the magnitude of the 
shift now also depends on n, the number of radial nodes. 
The origin of the dependence on sgn(m) may have a topo- 
logical interpretation. 



C. Inhomogeneous Magnetic Field 

We can also use a perturbative approach on eq. ([22]) 
to calculate the first order energy shift due to vector 
potentials with v > 1, which correspond to inhomo- 
geneous effective magnetic fields. The eigenf unctions 
and eigenvalues of the unperturbed Hermitian Hamil- 
tionian H (t)) = (l - p 2 ) V 2 - 2p^ are those of the 
harmonically trapped BEC when no effective magnetic 



field is present: Q(°> = p m P^ mm {l - 2p 2 )e im < i ' and 

= (2n 2 + 2n\m\ + 2n + |m|) 2 , where e' ^ is degen- 
erate with respect to the sign of m. The first order en- 
ergy corrections due to the perturbation H' — 7ml A| = 

arc t nc solutions of the secular equation 
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H' ++ - 4 
H'_ 



(i) 



H' 



H' _-e 



(i) 



with H', 



±± 



where + or — denotes 



the sign of m, and the constant be determined 

using Eqs. (fT5]) and (|2"3"]) . The off-diagonal terms vanish, 
allowing us to express the first order energy shift AE — 
E^ — E^ ' of the surface modes (n = 0) as 



AE = 



h 2 a £m 2{\m\ + 1) 



MR 2 (2| 



1)' 



(31) 



which is plotted in figure 3. The perturbative results 
agree well with our numerical calculations for the lowest 
energy excitations but breaks down as we increase n and 
m as would be expected since these excitations are no 
longer slowly varying in space. 

The observed shift on the excitation energy levels is 
reduced as we increase the exponent term v in the vec- 
tor potential from 1. This is not surprising as the radial 
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Frequency shift as a function of v {a.§t = 1) 



<] 0.05 





m=l (perturbative) 








m=2 (pertiubative) 








A oc r 1 






FIG. 3: Frequency shifts for the m = 1 and m — 2 modes 
due to the presence of the effective magnetic field as a func- 
tion of v, the radial exponent in the vector potential (Eq. 
©). The dotted lines correspond to results obtained solving 
Eqs. (|16|l - (|18| ) numerically and the solid lines to the results 
of the perturbative calculation (Eq. (I31|l . All frequencies are 
in units of the harmonic trap frequency uJt, and the parame- 
ters used were aol = —1 and 7 2 = 0.001, which describes a 
condensate well in the Thomas-Fermi regime. 



FIG. 4: Numerical results for surface mode (n = 0) excitation 
frequency as a function of the effective magnetic field strength 
ao£ for different vector potentials with radial exponent v. The 
lines correspond to the five lowest energy surface excitations, 
where the m — case is represented by the lowest frequency 
mode at ao£ = 0, and the m = 4 the highest. Only the 
v = 1 case exhibits energetic instabilities. Here ao = 0.5, 
7 2 = 0.001 and R = 1. The number of grid points used was 
128 and the system size was L=2.56 



coordinate p is < 1 and therefore for larger v the ef- 
fective magnetic field becomes more concentrated at the 
edge and permeates less to the centre of the cloud. As a 
consequence, the critical a§i for energetic instabilities to 
occur increases as the exponent v is increased. However, 
the perturbative treatment is valid only for << 1 
and so cannot be used to predict energetic instabilities. 
These instabilities can be inferred from the numerical cal- 
culations of the spectra for v — 2, 3, 4, presented in figure 
2. A full description of the numerical method is given in 
the next section. In contrast to the case of homogeneous 
effective magnetic fields we also note that the magnitude 
of the energy shift depends on the radial node index n. 
The effective magnetic field does not, however, affect the 
purely compressional (m — 0) modes with zero angular 
momentum. 



V. FULL EFFECTIVE TRAPPING POTENTIAL 
INDUCED BY THE LIGHT 

We now present the results of numerical calculations 
used to determine the excitations when the only assump- 
tion we make of our external trap is that it is harmonic, 
and include the additional effective trapping potential 
due to the interaction with the light and atoms. For re- 
alistic parameter values, with £ > 1 and ao < 1, the effect 
of the trapping potential induced by light eclipses that of 
the effective magnetic field, with profound implications 
for the excitation spectra. 

The exception is the case where v = 1. If ao is small, 
then the vector potential approximated by Eq. © cor- 
responds to a homogeneous magnetic field. The second 
term in Eq. (|10J) represents a uniform shift in the chem- 



ical potential throughout the cloud - an effect we can 
ignore by setting the effective trap minimum to zero - 
while the third term shifts the harmonic trap frequency 
downwards. The spectrum is then the same as Eq. (|26p 
except that we replace the trap potential ui t such that 

— Ms) 2 )'- <-> 



The spectrum exhibits energetic instabilities above a 
critical angular momentum of the light I, as in section 
IV. However, this is not the case for v > 1. In figure 4 we 
plot the dispersion curves obtained for v — 1, 2, 3, 4. An 
analytical approach is no longer possible when v > 1, cor- 
responding to an inhomogeneous effective magnetic field, 
and/or if the intensity ratio of the control and probe 
beam is such that the inequality a§ <C 1 is violated. 
In these cases we calculate the excitations solving Eqs. 
(fT5|) - (|T5]) numerically using the exact expressions for the 
potentials (Eqs. ([7][8])). First, we solve the radial version 
of Eq. (fig)) imposing first derivatives equal to zero at 
r = and r = L, with L the domain size. To do so we 
discretize the space and solve the set of coupled ordinary 
differential equations using a Newton method in which 
the derivatives are computed in Fourier space (20l . l2lj |. 
Solving Eqs. (fT6]) and (fl7|) then reduces to finding the 
eigenvalues of the discretized matrix associated to the 
linear problem (|16H17p . For v > 2 the mode frequencies 
only decrease to some value above zero until an interme- 
diate value of £, above which they increase monotonically 
with L 
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VI. SUMMARY AND DISCUSSION OF MAIN 
RESULTS 

In the studies of BECs subject to an external rotation - 
equivalent to applying a homogeneous effective magnetic 
field - the presence of instabilities in the excitation spec- 
tra has been shown to play a crucial role in determining 
the evolution of the condensate, particularly with respect 
to the nucleation of vortices. There are two types of in- 
stability which ought to be considered. Dynamic insta- 
bilities are associated with a complex excitation energy 
and departure from the initial configuration due to inter- 
action effects. The energetic instability relates to the ex- 
istence of excitations with negative energy, so that in the 
presence of dissipation the system can lower its energy by 
going into an 'anomolous' mode, and is a prerequisite for 
dynamical instability [22| . The dynamical instability has 
been credited as the primary mechanism for vortex nu- 
cleation in the rotating trap experiments [l7| . It has also 
been argued that while the dynamic instability helps in- 
duce vortex nucleation, the actual penetration of vortices 
in to the bulk is a consequence of the energetic instability 

The excitation spectra studied in this paper, both an- 
alytically and numerically do not exhibit any dynamic 
instabilities in that we do not observe a critical effective 
magnetic field at which the excitation energy becomes 
complex. Due to the increased complexity of adding an 
inhomogeneous effective magnetic field, however, we have 
restricted ourselves to consider trapping configurations 
which are radially symmetric. For the homogeneous ef- 
fective magnetic field dynamical instabilities are only ob- 
served with the introduction of an anisotropic trap [l7| , 
and this is likely to be the case for inhomogeneous mag- 
netic fields as well. 

Energetic instabilities, by contrast, occur readily. For 
all the vector potentials considered in section IV the nu- 
merical results display energetic instabilities beyond a 
critical field strength proportional to a^l. As we increase 
the radial exponent v in the vector potential we in turn 
must significantly increase the probe beam angular mo- 
mentum before the instability can be observed. When 
the full effective trapping potential is included the v = 1 
case exhibits an energetic instability but in the vast ma- 
jority of parameter space those for v > 1 do not. Only 
in the region where the maximum intensity ratio cxq ap- 
proaches 1 do we observe instabilities in the low-energy 
modes. However, it is in this region the adiabatic dynam- 
ics most easily breaks down [l[ and the model described 
in section II may become invalid. 

Another important feature of the spectra calculated in 
sections IV and V is the accidental degeneracies, which 
occur where the mode frequencies intersect in figures 1 



and 3. These degeneracies are likely to be manifested 
through transfer of excitations from higher order modes 
to lower order modes at the point of degeneracy. This 
kind of phenomenon could have implications for the con- 
densate evolution depending on whether the effective 
magnetic field is switched on instantly or its strength 
is adiabatically ramped up to a final value. For ex- 
ample, an anisotropically trapped condensate in a ho- 
mogeneous effective magnetic field naturally undergoes 
m = 2 quadrupole oscillations, but the higher order 
modes exhibit dynamic instability [l7| . eventually lead- 
ing to vortex nucleation. To properly account for the 
mode-coupling we would need to move beyond Bogoli- 
ubov - de Gennes theory which does not account for in- 
teractions between the degenerate modes [TBI ]. 

With inclusion of the full effective trap induced by the 
interaction with the light we reach a scenario where there 
are no accidental degeneracies as v is increased beyond 
1. The vector potential appearing in the Bogoliubov-de 
Gennes equations is less important in determining the 
excitation frequencies in comparison with the effective 
trap. A direct measure of the contribution of the vector 
potential is the splitting of the + and - m modes. For the 
potentials with v > 2 studied in section V, we find that 
the splitting is typically < 5% of the actual mode fre- 
quency if no effective magnetic field were present. This 
is in stark contrast to section IV where the splitting can 
easily exceed twice the mode frequency causing an ener- 
getic instability for realistic parameter values. 

In this paper we have deliberately paid particular at- 
tention to the surface modes (n = 0, m ^ 0), which are 
naturally excited by adding the effective magnetic field to 
an anisotropically trapped condensate (23|. We predict 
that for inhomogeneous effective magnetic fields vortex 
nucleation will be inhibited due to the dominance of the 
trapping potential induced by the light. In general, it is 
the ratio of the cyclotron frequency to the effective trap- 
ping frequency which drives the dynamics of the system 
and both of these depend on a^l. It is therefore often 
desirable to introduce a counter potential to act against 
the additional trapping terms due to the interaction with 
the light, as described in section IV. 
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